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ABSTRACT 



We present new JHK spectroscopy {R ~ 5000) of GQ Lup b, acquired with 
the near-infrared integral field spectrograph NIFS and the adaptive optics system 
ALTAIR at the Gemini North telescope. Angular differential imaging was used 
in the J and H bands to suppress the speckle noise from GQ Lup A; we show that 
this approach can provide improvements in signal-to-noise ratio (S/N) by a factor 
of 2 — 6 for companions located at subarcsecond separations. Based on high- 
quality observations and GAIA synthetic spectra, we estimate the companion 
effective temperature to T^q = 2400 ± 100 K, its gravity to log (7 = 4.0 ± 0.5, and 
its luminosity to log(L/L0) — — 2.47± 0.28. Comparisons with the predictions 
of the DUSTY evolutionary tracks allow us to constrain the mass of GQ Lup 
b to 8 — 60 Mjup, most likely in the brown dwarf regime. Compared with the 
spectra published by Seifahrt and collaborators, our spectra of GQ Lup b are 
significantly redder (by 15 — 50%) and do not show important Pa(3 emission. 
Our spectra are in excellent agreement with the lower S/N spectra previously 
published by McElwain and collaborators. 

Subject headings: planetary systems - stars: low mass, brown dwarfs - stars: pre-main 
sequence - stars: individual (GQ Lup) - techniques: spectroscopic - techniques: high 
angular resolution 
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Introduction 



An important goal of exoplanet imaging is understanding how planetary systems form 
and this requires observing planets when they are very young, i.e. at ages less than a few 
million years. Ten years after the first discovery of an extra-solar giant planet (EGP), the 
imaging era of exoplanet science has finally started. The year 2008 has been particularly 
prolific in tha t matter with the d iscovery of the first planetary system around the ~ 60 Myr 



old HR 8799 flMarois et al. 



Myr old Fomalhaut (Kalas et al 



( Lagrange et al 



( iLafreniere et al. 



2008h . of th e exoplanet candidates in orbit around the 100 — 300 



20081 ). in orbit around the ~ 12 Myr old j3 Pictoris 



20081) . the potential companion to the ~ 5 Myr old IRXS J160929-2 10524 



20081 ) and the planetar y mass obiect 



0.9 to 3 Myr old T-Tauri star CT CHA ( 



PMO ) candidate evolving round the 



Schmidt et al. 



20081 ). These new discoveries add to 



a few previous detections of planet candidates suc h as the PMO in orb it around the young 



10 Myrs) 25 Mj^p brown dwarf (BD) 2M1207 fIChauvin et al. 



2004r) and the com panion 



in orbit around the ~ 1 — 10 Myr old T Tauri star GQ Lup (iNeuhauser et al. 



20051). 



GQ Lup b is particularly interesting; the co-moving companion is relatively bright 
(K=13) and lies at 0.7 arcsec (100 ± 50 AU) from the 7th magnitude K7eV primary. 



Despi te several studies, the mass of GQ Lup b is still poorly constrained. 



Neuhauser et al. 



( I2OO5I ) identified prominent CO and H2O bands typical of cool M9— L4 dwarfs but could 
not constrain accurately the mass of the companion (1 — 42 Mjup) because the effective 
temperature (Tpff ; 1600 — 2500 K), t he gravity and the distance of the object were not 



well determined. 



Marois et al. 



m 



fl2007f ) have analyzed archived HST WFPC2 and NICMOS 



data and Subaru CH4, H, Kg and L' band images of GQ Lup b and provided estimates 
of its bolometric luminosity (log(L/L0) = —2.42 ± 0.07 at 140 pc), its radius (0.38 ± 
0.05 Rq) and its effective tem perature (2335 ± 100 K). Assuming a log g = ?> and based 



on the evolutionary models of 



Baraffe et al. 



(120031 ). they inferred a mass ranging from 
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10 to 20 Mjup. iMcElwain et al.l (120071 ) secured J- and if-band spectra and found the 
companion to be somewhat hotter than previously reported with a spectral type between 
M6 and LO corresponding to effective temperatures between 2400 and 2900 K assum ing 
the effective temperature - spectral type relationship fro m 



Using the evolutionary models from 



Burrows et al. 



Golimowski et al 



fll997l ) and 



Chabrier et al. 



(2004). 



(120001) and 



their derived bolometric luminosity of log(L/L0) = —2.46 ± 0.14 they infer r ed a mass 



of 10 — 40 Mjup. A mass estimate has also been reported by 



Seifahrt et al 



mm who 



acquired H- and /T-band sp e ctra a nd compared them to the predictions of the GAIA 



model from 



Brott fc Hauschildt 



(120051 ). The best fit resulted in Teg = 2550 - 27507^ and 
log g = 3.7 ± 0.5 dex. Fror n these values, their estimated log(L/Lo) = —2.25 ± 0.24 and 



by comparison with 2M0535 (IStassun et al. 



20061 ) which has a dynamically determined mass 



of 36 ± 3 Mjup and a similar age as GQ Lup, they derived a mass between 4 and 36 Mjup. 
A summary of previously reported physical parameters of GQ Lup b is given in Table HI 

In light of the International Astronomical Union (lAU) working definition which states 
that a planet must have a mass below the deuterium burning limit of ~ 13 Mjup, the 
nature of GQ Lup b thus remains uncertain. Further progress will be made through a 
better distance estimation of the system and better spectroscopic characterization of the 
companion. 

In this paper, we present new R ~ 5000 JHK spectroscopic observations of the 
GQ Lup system taken with the integral field spectrograph NIFS in combination with the 
adaptive optics system ALTAIR at the Gemini North telescope using angular differential 
imaging (ADI) to improve the signal-to-noise ratio (S/N) of the companion. The resulting 
spectrum is analyzed and compared with synthetic atmospheric spectra and evolutionary 
tracks to infer the physical parameters of the companion. 
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Observations 



The GQ Lup system was observed in tlie J, H and K bands respectively on the nights of 
2007 May 30, June 1 and June 27 at the Gemini North Telescope (p rogram GN-2007A- O-46) 



using NIFS combined with the ALTAIR adaptive optics system (IHerriot et al. 



19981) 



NIFS is a near-infrared integral field spectrograph (IFS) based on the image sheer concept 



(IMcGregor et al. 



2002I ). It has a spectral resolving power of 5000 in each band, a total field 



of view of 3.0" x 3.0" and rectangular spatial resolution elements of 0.04" x 0.1". 

All observations were obtained with the Cassegrain rotator turned off to allow rotation 
of the field around the host star. The J- and if-band sequences consisted of 18 x 400s 



and 27 x 300s exposures, respectively, without an y dithering 



Marois et al. 



his a. 



20061 ) 



lo wed enough fielc 



Lafreniere et al. 



rotati on to use angular differential imaging (ADI; 
(120071 )) to suppress the speckle noise. As the PSF is better sampled and AO achieves better 
performance in the K band, observations were obtained with a 9-position square dithering 
pattern with each position offsetted by 0.4"; in total 72 exposures of 50s were obtained. 
The median seeing was measured to be 0.44", 0.61" and 0.67" in J, H and K respectively 
and the observations were conducted under a median air mass of 1.79 in both J and H and 
of 1.83 in K. A spectroscopic standard for telluric correction was observed at five random 
dither positions having a maximal excursion of ±0.7" after the J-band sequence and both 
before and after the H and K bands observations. The spectroscopic standards for telluric 
correction were HIP 82714 in J and HIP 69021 in both H and K. Sky and arc lamp frames 
were taken at the end of each sequence. Darks and flats as well as Ronchi flats allowing the 
calibration of the spatial distortion of the camera were taken during daytime calibration. 
No astrometric calibration data (plate scale and orientation) were obtained. An observation 
log is shown in Table [2l 
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3. Data Reduction 

3.1. Primary data reduction and data cube reconstruction 

The data cube reconstruction was done with the IRAF Gemini data reduction package 
version 1.9.1a. The calibration data were first processed to construct a fiat field, a bad pixel 
map, a wavelength solution and a map of the curved spectral traces. 

The telluric data cubes were then reconstructed. The images were rectified to correct 
for the trace curvature and re-sampled to yield a common constant wavelength step for 
all spatial elements. A sky cube was built from the median of the dithered spectroscopic 
standard data and subtracted from each telluric cube. The latter were then corrected for 
atmospheric dispersion by fitting a Gaussian kernel to the PSF to find its center in each 
spectral slice and by translating the PSFs to the center of the field-of-view. The median 
telluric spectrum was computed and divided by the spectrum of a blackbody at 10 000 K, 
the approximate AOV spectroscopic standard temperature. 

The final step consisted in reducing the science data and building their associated 
data cubes. The sky frame obtained after the science observations was first subtracted 
from each image of the time-lapsed sequence. The science data cubes were extracted 
following the same procedure used for the telluric data. Telluric correction was then 
applied. Atmospheric dispersion was then corrected as before. Interstellar extinction was 
considered to be ne gligible as most of the extinction associated with GQ Lup A is believed 



to be circumstellar ( iMcElwain et al. 



20071 ). The rectangular IFS spatial resolution elements 
were interpolated to yield a spaxel size of 0.043" x 0.043" in the final reconstructed cube. 
The term spaxel will be used to refer to a single {x, y, A) element in the final data cube. 
This is applied to all IFS images to get a final dataset composed of time-lapsed data cubes. 

An example of a collapsed data cube, i.e. integrated over the full wavelength range 
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of the photometric band, obtained for each observed band is shown in the top panels of 
Figure [TJ A bright vertical stripe aligned with NIFS slitlets is present in all our data. NIFS 
is known to have some scattered light in each independent channel which can result in an 
enhanced sky proportional to the star brightness in each slitlet and could explain such a 
bright stripe. 

In preparation for the extraction of the signal of GQ Lup b, the bright signal from the 
primary PSF was subtracted to first order by applying an unsharp mask with a Gaussian 
kernel of six pixels FWHM, followed by the subtraction of a median radial intensity 
profile about the center of the PSF. The images were then convolved by a 1 pixel FWHM 
Gaussian kernel to attenuate the pixel-to-pixel noise resulting from the interpolation of the 
rectangular IFS spatial resolution elements to the square data cube spaxels. The resulting 
residual broad band images obtained after collapsing the data cubes are shown in the lower 
panels of Figure [H Note that the residual vertical stripe is subtracted by the ADI speckle 
suppression algorithm (see next section) and does not impact the final J and H spectra. 
In the K band, however, no such algorithm is used. We therefore only considered the first 
45 data cubes obtained, in which the companion does not overlap with the vertical stripe. 
The K band spectrum was extracted from these residual images by summing the fiux in a 
2 spaxel diameter circular aperture. 



3.2. Angular Differential Imaging 



As can be seen in Figure [T], the residual speckle noise is significant at the separation of 
the companion and this severely limits the accuracy by which the spectrum can be extracted. 
In principle, thi s speckle noise could be attenuated using spectra l deconvolution tech niques 



(ISoarks k FordI 12002: iMcElwain et all 120071 : iThatte et all 120071 : iJanson et al.l l2008l : 
Lavigne et al.l 12009. in preparation! ) but, in practice, the companion radial displacement 
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is only 2 — 3 NIFS spaxels in a data cube rescaled such that the speckles at different 
wavelengths overlap. This, combined with the companion brightness, makes an accurate fit 
of the speckle pattern contaminating the companion spectrum very hard to obtain. 



Marois et al. 



he approach chosen instead was to use the ADI technique developed by 
( I2OO6I ). This technique consists in acquiring a series of images with the Cassegrain rotator 
of an alt-az telescope turned off. This allows field rotation between each observation of the 
sequence. The companion is then rotating around its host star in the time-lapse sequence 
while the quasi-static speckles introduced by the defects of the telescope optical elements 
stay fixed. For each data cube of the time series, a reference data cube can be built from 
the other data cubes, and used to subtract the speckle pattern and attenuate the speckle 
noise. 

Examples of images obtained after the subtraction of their median radial profile taken 
at intervals of 65 minutes in the H band are shown in Figure [2l This corresponds to an 
image at the beginning, the middle and the end of the whole sequence. The long speckle 
life time is clearly seen in that sequence while the companion revolves around its host star 
due to field rotation. 



The reference data cube wa s built using the "loca. 



(LOCI) algorithm developed by 



Lafreniere et al 



ly optimized combination of images" 



(I2OO7I ). This algorithm builds a reference 



PSF having the highest speckle correlation with respect to the analyzed frame from the 
dataset composed of the time sequence of images. This is done by first dividing the frame 
in optimization subsections. The reference PSF is then built for each zone by computing 
the linear combination of images that guarantees the lowest possible residual noise after 
subtraction. Only the images of the time sequence in which the companion had a sufficient 
displacement compared to the analyzed frame are used to build the reference PSF to avoid a 
significant companion flux loss. The resulting minimal time separation between two frames 
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is the one yielding the best companion S/N when two frames are differentiated. This was 
found to be 2400s in J and 2700s in the H band which respectively corresponds to 6 and 9 
images in the time sequence and to a minimal companion displacement of 2.7 and 2.9 pixels. 

Since the companion position is already known, a single horseshoe shaped optimization 
zone centered on the host star and having a mean radius corresponding to that of the 
companion was used. The radial width is taken to be 4 pixels. The missing azimuthal zone 
is the one in which the companion is moving through the time sequence. This algorithm is 
applied to each spectral slice of each data cube. The resulting data cubes are then rotated 
to align the field of view and summed to build the final ADI data cube. The companion 
spectrum was extracted from this final cube by summing the flux in a 2 spaxel diameter 
circular aperture. 

Applying this ADI process to each data cube has the advantage of increasing the 
S/N of the companion spectrum by significantly attenuating the speckle noise. However, 
a fraction of the comp anion signal is lost through the application of the LOCI algorithm 



( iLafreniere et al 



20071 ). Moreover, the fact that each spectral slice is treated independently 
from the others can have an impact on the resulting spectral shape. Hence, we have 
developed an algorithm to determine the influence of ADI on the companion spectrum. 

For simplicity, let us consider only a single slice of a given data cube in the sequence, 
keeping in mind that the procedure described below is repeated for each slice and each cube. 
Following the execution of the above LOCI subtraction algorithm, the residual integrated 
flux at the companion position in the cube considered can be written as 



F = I-Y,C^k, (1) 



where / is the integrated flux in the slice considered before LOCI subtraction, Jj is the 
same but for cube and {cj} is the set of coefficients found by the LOCI algorithm. The 
integrated flux / can be expressed as the sum of a signal C coming from the companion and 
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a signal S coming from quasi-static speckles, yielding 

F={S + C)-J2c^iSi + Q) = S,^ + C-J2 CiC^, (2) 

i i 

where S'res = "S* — Ylii (^i^i is the residual speckle signal, which on average will be close to 

zero. Our goal is to find the true companion signal C, and so we must estimate the quantity 

Y2^CiCi, i.e. the fraction of the companion signal that was subtracted by the algorithm. 

This was done by introducing a fake companion of known fiux and fiat spectrum into the 

images, at the same separation as GQ Lup b but 180° from it. The fake companion was 

taken to be a shifted and intensity scaled-down version of the primary star; its PSF should 

thus be very similar to that of GQ Lup b. The LOCI algorithm was applied to the sequence 

of data with the fake companion by using the exact same coefficients as for the data without 

the fake companion. Then the fake companion residual fiux F' was retrieved. Since the 

true fake companion signal C and the residual speckle signal S^^^ at the location of the 

fake companion are known (the latter from the execution of the algorithm without the fake 

companion), the fraction of the fake companion fiux that was subtracted by the algorithm 

in a given spectral slice can be readily computed as 

_ Si CjQ' S'j.^^ + C' — F' 
JA = —Qi— • y-^) 

This fraction should be the same for GQ Lup b as for the fake companion since the 
same coefficients were used for both and they have the same PSF. Note also that f\ is 
independent of the fake companion azimuthal position as ^ • CiC[ only depends on the 
relative angular position of the companion PSF between the analyzed frame and the cube 
i. So, the true signal of GQ Lup b can be calculated as 

C = (4) 

where in practice S'res is neglected; the resulting error on C is computed by measuring the 
residual speckle noise at the companion radius. The correction factors f\ computed for the 
J and H ADI spectra by the above procedure are shown in Figure [3l 
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Finally, to make sure that the above procedure did not introduce any artificial feature 
in the spectra, mainly in their continuum shape, we have performed a simple spectrum 
extraction directly from the reduced data cubes without any ADI processing. The resulting 
spectra were compared with the spectra obtained with ADI processing and no significant 
difference was found, save for a lower S/N. More precisely, a similar procedure as for the 
K band data was used to extract a non-ADI spectrum from the data cubes in which 
the companion is the least affected by the vertical stripe in the J and H bands. The 
resulting ADI spectrum was then divided by the non-ADI one to characterize the continuum 
dispersion. The median non-ADI J band spectrum of the first six data cubes had only a 
5% continuum difference with the ADI spectrum. Similarly, the first eight and last four 
data cubes in the if-band were used to compute the median non-ADI spectrum and the 
resulting continuum had only a 6% slope difference with the ADI spectrum. 

4. Results 

4.1. Speckle noise attenuation performance 

The ability of the algorithms described in the previous section to attenuate the speckle 
noise in NIFS data is evaluated in this section. The speckle noise attenuation factor is 
defined as the ratio of the noise in the initial image over the noise in an ADI image. The 
noise in ADI images was corrected for the companion fiux loss during the ADI process using 
fake companion implants. Each spectral slice was evaluated independently. Wc report the 
median value of speckle noise attenuation in all the spectral slices as a function of angular 
separation from GQ Lup. The noise was evaluated by computing the standard deviation of 
spaxels within an annulus of one spaxel width at each radius. GQ Lup b was masked out 
from the images. 
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The results are plotted at the top of Figure H] for the J and H bands while the contrast 
curves achieved are plotted at the bottom of the figure. For a single data cube, the ADI 
subtraction provided a noise attenuation by a factor of 1.5 to 2 in both J and if, while 
the combination of all ADI data cubes after alignment of their field of view provided an 
additional noise attenuation by a factor of 2 to 3. Overall, the speckle noise is attenuated 
by a factor of 3 — 6 in J and a factor 2 — 5 in iJ. 



4.2. Photometry 



The relative photometry of the companion compared with its host star is first evaluated. 
Unfortunately, the peak of GQ Lup A is saturated or in the detector non-linear regime 
over most of the J band and the entire H band, preventing us from calculating differential 
photometry over those bands. In J, over the narrow spectral range in which GQ Lup A is 
in the detector linear regime (1.31 — 1.35 /im), we calculated a contrast of 6.53 ± 0.06 mag 
for the companion. 

The uncertainty was measured by computing the standard deviation of the flux in 20 
apertures identical to the one used to integrate the companion flux at the same radius as 
the companion. The peak of GQ Lup A is below non-linearity in the entire K band; we 
measured a differential photometry of AKg = 6.35 ± 0.11 mag. 



GQ Lup A is a known variable classical T Tauri star (e.g. 



Batalha et al. 



Broeg et al. 



(1200 ih I 



(120071 ) reported it to vary by ±0.44 r nag and ±0.22 rn a g in J and Kg bands 



respectively over a period of 8.38 ± 0.2 days while 



Neuhauser et al. 



(120081 ) found it to vary 



by a maximal amplitude of ±0.20 mag in the K,, band between 



Using the GQ Lup A photometry reported by iNeuhauser et al. 



Mav_ 2005 and Feb 2007. 



(l2008f ) of K, = 7.18 ± 0.20 



mag and the J-band photometry from 2MASS J = 8.605 ± 0.44 with the updated error 
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bar from 



Broeg et al 



(120071 ). we find a GQ Lup b photometry of Ks = 1 3.45 ± 0.25 and 



J = 15.13 ± 0.44. The value is consistent with = 13.39 ± 0.08 from 



( 2008 ) and the J-band value is consistent with 14.90 ± 0.11 reported by 
J2007f ). 



Neuhauser et al 



McElwain et al. 



4.3. Spectroscopy 



Our spectra of GQ Lup b are shown in Fig. 5. The S/N are 145, 76 and 15 at 
R ~ 5000 in JHK respectively where the noise was computed from the dispersion of the 
fiux measured in 20 apertures situated at the same radius as the companion. The integrated 
fiux in each aperture were corrected by f\ (see section [321) • 

In the same figure, our J HK spectrum of GQ Lup b is compared with th e previously 



published results obtained by 



Seifahrt et al. 



mm and 



McElwain et al 



noted feature is that our spectra are redder than the ones reported by 



2007 



The first 



Seifahrt et al 



(l2007f ) 



in all three ban d s whi le they have a similar continuum slope as the ones reported 



McElwain et al. 



(120071 ). Dividing our spectra by the ones reported by 



Seifahrt et al. 



by 



mm 



reveals differences in the overall spectral slopes of 22%, 48% and 15% res pectively in J, H 



and K . This difference is of 8% in J and of 17% in H when compared to 



McElwain et al. 



( 120071 ) where most of the if-band difference lies in the bluer spectral range with a 14% 
continuum difference below 1.63 yum. The difference between the two spectra at the redder 
end is then of only ~ 3%. 

A second difference is seen for the Pa(3 emission line, which is weaker by about an 



^Note that the 



McElwain et al. 



(120071 ) spectra were taken during the OSIRIS commis- 
sioning run. The lower S/N is not representative of the camera performance as several 
adjustments were made to the camera after these observations were taken. 
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order of magnitude in our spectrum compared to that of ISeifahrt et al.l (120071 ). While it 
is detected with an equivalent width of EW(Pa/3) = —3.83 ± 0.15 by these authors, it 
is detected with an EW{Paj3) = —0.46 ± 0.08 in ours. Emission-lines of PMS accreting 
objects are known to vary either due to accretion variability, to projection effects on the 
accreting column rotating around the object or to hot spots that are not a lways visible to 



the observer due to the object rotation (see iScholz fc Jayawardhanal (120061 ) and reference 



th erein for more detai ls). The strongest PajS variations found in the literatur e were repor ted 



by 



Natta et al 



Gatti et al 



(120061 ) who analyzed a sample of 14 BDs also observed by 
(l2006l ) one to two years later. They found a variation of a factor of two at most for all the 
objects in their sample except for one that varied by a factor of three. This makes GQ Lup 
b significantly more variable than any other BD previously observed. 



This fac t along with t 



spectra from 



l e con tinuum differences raise the possibility that GQ Lup b 



Seifahrt et al 



(120071 ) are contaminated by some host star flux. Indeed, such a 
contamination could expla in both the bluer con tinuum and the signiflcantly stronger Pa(3 
emission line measured by 



Seifahrt et al. 



(120071 ) since the spectrum of A is bluer than B 



and exhibits strong Pa/3 emission. However, the similarities between the KI doublets EW 
reported in Table [3] for both spectra argues against this. 

The use of PSF subtraction algorithms could potentially introduce artefacts in the 
companion spectrum. The fact that different algorithms were used in the three compared 
papers could then explain the discrepancy between their spectra. As discussed before, we do 
not believe that artefacts were introduced by our ADI subtraction procedure given the small 
difference (respectively, 5 and 6% for the J and the H band) between the median non-ADI 
and ADI spectra. This assertion is rei nforced by the prese nce of the same discrepancy 



between our results and the ones from ISeifahrt et al. 



(120071 ) in the K band where no ADI 



processing was done. At this point, we do not have any satisfactory explanation for the 
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discrepancy between the results of the different papers. 



5. Discussion 



In Figure [6], our spectra are compared with the ~ 5 Myr late-M to ear ly-L dwarfs 



USco J160830-233511, USco J161302-212428 and USco J163919-253409 from 



the 1 - 50 Myr 2MASS J014158-463357 (IKirkpatrick et aL 



LodieuetaL 



20061 ) and some 



older field dwa.rfs ta ken from the NIRSPEC Brown D warf Spectroscopic Survey (BDSS) 



(IMcLean et al. 



2003 



and from 



Gushing et al. 



(2005 



A first feature that is seen by 



comparing GQ Lup b with the field dwarfs is that the absorption from H oO at 1.34 p,m 



Gorlova et al. 



known as being mostly temperature dependent in low-mass objects (e.g. 
(120031 ) ). is similar to objects of spectral types between LO and L2. This is also supported by 
the comparison of the continuum slope shortward of ~ 1.55 fim and longward of ~ 1.68 fim, 
while the region in between has been reported as being particularly affected by the reduced 
H9 coUisionally- induced absorption (GIA) in young, low gravity objects such as GQ Lup b 



teorysow et al. 



1997 



Kirkpatrick et al. 



2006|). This also explains the better agreement in 



that zone with the younger objects. H2 GIA has its strongest effect in the K band where 
GQ Lup b does not exhibit the decrease in flux seen between 2.15 and 2.30 fim in field 
objects and has significantly shallower GO absorption. The i^'-band spectrum is better 
reproduced by the young objects, with the LI spectra being the most similar. 

Hence, we find a spectral type of Llifcl for this object. This is at the cooler end of 



the M6-L0 spectral type reported by 



McElwain et al. 



(I2OO7I ) and is consistent with the 



^The spectra can be found at http:/www. iac.es/galeria/nlodieu/ 

^The spectra can be found at http:/www. astro. ucla.edu/~mclean/BDSSarchive/ 

^The spectra can be found at |http : / / irtfweb . if a. hawaii . edu/ ~spex/ spexlibr ary / IRTFlibr ary. html 
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M9-L4 spectral type inferred bv lNeuhauser et al.l (120051) . This is verified by computing tfie 



H2O index at 1.5 fim defined by 



AUers et al. 



(l2007l ). which is beheved to be independent 
of surface gravity. We find a spectral type of LO ±1 also consistent wit h our spectral 



type estimate. The Tes - spectral type relationship of 



Golimowski et al. 



(12004h implies 



Teff = 2250 ± 134 K. However, this relation was only calib rated for evolved f i eld o 



and thus should be used with caution with young objects. iKirkpatrick et al. 



Djects 



fl2008h also 



demonstrate that objects of a similar optically determined spectral type exhibit a significant 
dispersion in their near-IR relative color. This illustrates the challenge of using such a 
relation to infer Tes, especially from near-IR data. 

One can also estimate the effective temperature of GQ Lup b and its gravity by 
comparing our JHK spec tra with synthetic ones generated by the GAIA model v. 2. 6.1 



(iBrott fc Hauschildt 



20051 ). The grid explored includes R=5000 spectra ranging from Teg 
= 2000 to 2900 K and from log g = 1.5 to 5.5 dex with increments of 100 K and 0.5 dex 
respectively. A minimization a lgorit hm was used to find the best fitting parameters. 



As reported by iKirkpatrick et al. 



(120061 ). current atmosphere models poorly reproduce the 
strong H2O absorptions between 1.32 and 1.60 /im and between 1.75 and 2.20 /im. Hence, 
the H band was not used in the fit and those zones were masked out of the J- and i^-band 
spectra. The Pa/3 emission line in the J band was also masked out of the fit. The fT-band 
fit is poor and leads to an optimal Tgff > 2900 K which is largely inconsistent with any 
previous T^^ estimates for this object. However, our i^'-band spectrum is very similar to 
other young brown dwarfs of similar spectral types in the zone fitted as seen in Figure El 
This indicates that there is still some discrepancy between models and young low-gravity 
objects in that band. Hence, only the J-band spectrum was used to estimate Teg and log g. 
Good fits were obtained for Tcfj = 2300 — 2500 K and log g = 3.5 — 4.5 with a best fit 
value of Teff = 2400 K and log g = 4.0. The model and observed spectra are compared in 
Figure [7| for the J band. The observed H- and i^-band spectra as well as the synthetic 
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enes obtained for those best fitting values are shown in Figure [H Hence, we adopt Teg = 
2400 ± 100 K and log g = 4.0 ± 0.5 dex. The same exercise was repeated with grids of 
AMES-dusty and AMES-cond synthetic spectra and all converged to the same solution. 

From th e se res ults, the mass of GQ Lup b is evaluated following the procedure used in 



Marois et al. 



(120081 ). The model spectra fiux is first adjusted to best fit the photometric 
data points reported in this work and in the literature through a minimization algorithm. 
The resulting spectra along with the photometric measurements used for the fit are plotted 
in Figure [9l The luminosity is then evaluated to be equal to log(L/L0) = —2.47 ± 0.28 
dex by integrating the model fiux over the whole wavelen gth range and by using the 



companion distance es timate from 



Neuhauser et al. 



McElwain et al 



fl2008h of 139 ± 45 pc. Adopting the 



( 120071 ) age estimate of 1 — 10 Myr f or the companion and plotting the new 



luminosity on the DUSTY evolutionary tracks from 



Chabrier et al 




(1200 



yields a mass 



estimate between ~ 11 — 60 Mjup (see Figure [TO]) putting GQ Lup b most likely in the 
brown dwarf regime. 

Even though e v olutionary mo d els are unc ertain at the young age of GQ Lup 



( iBaraffe et al. 



2002) 



Stassun et al. 



(12006 



20071 ) provide a first anchor point with the 1 



+2 
1 



Myr eclipsing binary 2M0535-05. The lower mass object of this system (36 ± 3 Mjup) is 
particularly interesting as it has a mass comparable to the one estimated for GQ Lup b. 
The comparison of its dynamically measured mass wit h the prediction of the evolutionary 



Chabrier et al 



(I2OOOI ) evolutionary 



model gives an estimate of the model systematic error, 
tracks predict a luminosity of log(L/LQ) = —2.00 for an object with such a mass and then 
underestimate it by a factor of ~ 1.5. Even though this lies only at ~ 1.5a of its measured 
luminosity and, hence, does not constitute a convincing evidence of a model systematic 
error, it is still interesting to include this error in our analysis to take into account the 



^The tracks were retrieved from 



http: //perso. ens-lyon.fr/isabelle.baraffe/ 
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model uncertainty. Updating our mass estimate with this systematic error leads to a revised 
lower mass value of ~ 8 Mjup which is below the boundary between a planetary mass object 
and a brown dwarf. Our final mass estimate is then of 8 — 60 Mjup for GQ Lup b. 

6. Conclusions 

In this paper, we first derived a method to characterize and compensate for the effect of 
the ADI speckle suppression algorithm apphed to IFS data. We found that the ADI process 
provided a speckle noise attenuation by a factor greater than 3 in the J band and greater 
than 2 in the H band. We extracted spectra in the J, H and K bands for the low-mass 
companion GQ Lup b. The companion Pa/3 emission line was marginally detected. The 
extracted spectra allowed us to constrain its spectral type to LI ± 1 and to re-evaluate its 
Tgff to 2400 ± 100 K and log g to 3.5 — 4.5. This leads to the re-evaluation of the bolometric 
luminosity to log{Li,oi/ Lq) — —2.47 ± 0.28 and to a mass estimate of 8 — 60 Mj^p by 
comparison to predictions from evolutionary models. This defines GQ Lup b as being most 
likely a brown dwarf. 

Based on observations obtained at the Gemini Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc, under a cooperative agreement 
with the NSF on behalf of the Gemini partnership: the National Science Foundation 
(United States), the Science and Technology Facilities Council (United Kingdom), the 
National Research Council (Canada), CONICYT (Chile), the Australian Research Council 
(Austraha), Ministerio da Ciencia e Tecnologia (Brazil) and SECYT (Argentina). 

We would like to thank Andreas Seifahrt, Ralph Neuhauser, Michael McElwain, James 
Larkin and J. Davy Kirkpatrick for sharing their spectra with us. We would also like to 
acknowledge Eric Steinbring at the Canadian Gemini Office for his support during the time 
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Table 1. GQ Lup b previously published results 



Neuhauser et al. (2005) Marois et al. (2007) McElwain et al. (2007) Seifahrt et al. (2007) This work 



Age (Myr): 


1 ± 1 


1 ± 1 


1-10 


1±1 


1-10 


Distance (pc): 


140 ± 50 


140 ± 50 


150 ± 20 


140 ± 50 


139 ± 45 


Teff (K): 


2050 ± 450 


2335 ± 100 


2600 ± 300 


2650 ± 100 


2400 ± 100 


log g: 


2.52 ± 0.77 


3 (assumed) 




3.7±0.5 


4.0 ±0.5 


Spectral type: 


M9 - L4 




M6 - LO 




Ll±l 


Radius (-Rjup): 


~ 2 


3.7 ± 0.5 








log(L/L0): 


-2.37 ±0.41 


-2.42 ±0.07 


-2.46 ±0.14 


-2.25 ± 0.24 


-2.47 ±0.28 


Mass (A^jup): 


1-42 


10 - 20 


10 - 40 


4-36 


8-60 



Table 2. Observation log 



Band 


Date 


Number of 


Exposure 


Median 


Median 


Total field 


Dither 






exposures 


time (s) 


seeing (") 


airmass 


rotation (°) 


(") 


J band: 


2007 May 30 


18 


400 


0.44 


1.79 


34.3 


No 


H band: 


2007 Jun 1 


27 


300 


0.61 


1.79 


38.8 


No 


K band: 


2007 Jun 27 


72 


50 


0.67 


1.83 


33.1 


±0.7 



Table 3. Emission/absorption line equivalent width comparison 



Source 


EW[Pa/3] 


EW[KI(1.169 fim)] 


EW[KI(1.178 fim)] 


EW[KI(1.244 fim)] 


EW[KI(1.253 fim)] 




[A] 


[A] 


[A] 


[A] 


[A] 


This work: 


-0.46 ± 0.08 


2.12 ± 0.32 


2.97 ±0.44 


2.30 ±0.26 


1.40 ±0.17 


Seifahrt et al. (2007): 


-3.83 ±0.12 


2.22 ±0.19 


2.51 ±0.40 


1.72 ± 0.15 


1.61 ±0.17 
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Fig. 1. — Detection of the companion to GQ Lup in J, H and K from left to right. Each 
image has a field of view of 3.0" x 3.0". The top, center and bottom panels respectively show 
examples of collapsed data cubes before and after the subtraction of a radial profile and 
after the application of the ADI algorithm. Arrows indicate the companion position. Each 
image was divided by the host star flux integrated in a two pixels diameter circular aperture 
and the same dynamic range is used to display all bands. Images after the subtraction of 
the radial profile and after the application of the ADI algorithm are also shown on the same 
dynamic range. 




Fig. 2. — if-band time sequence of images with the Cass-rotator turned off. Images are 
separated by 65 minutes each. The speckles long life time is clearly seen from this sequence. 
Arrows indicate the companion position. The field has rotated by 21.0° in the second image 
and by 37.7° in the third. 
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1.20 1.25 1.30 

Wavelength (|im) 




1 .50 1 .55 



1.60 1.65 1.70 
Wavelength (|im) 



1.75 



Fig. 3. — Correction factors applied to the extracted J- and if-band spectra at a spectral 
resolution of i? ~ 5000 to compensate for the flux loss during spatial filtering and the ADI 
process. The fake companions used to compute this factor had a flat spectrum (see section 
13.21 for more details). 
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Fig. 4. — The speckle noise attenuation performance is plotted in the upper panels for 
collapsed data cubes in J and H bands. The dashed Hnc represents the improvement between 
a single image of the time sequence and its associated ADI image. The solid line is the 
attenuation obtained between a single image and the final ADI image constructed by aligning 
the field of view and adding all the single ADI images. The bottom panels show the median 
detection limit in the non-ADI images (dot-dashed line), the median detection limit in single 
ADI images (dashed line) and after field of view alignment and addition of all the single ADI 
images (solid line). The contrast or detection hmit is defined as the Scr speckle noise level 
in an image. The ADI curves were corrected for the companion flux loss resulting from the 
reference image subtraction. Note that the reference image was constructed with a singular 
LOCI optimization zone situated at the companion radius and, hence, the reported curves 
are the lowest achievable performance with NIFS at other radii. 
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Fig. 5. — Comparison of GQ 
ones previously published by 



up b JHK spectra obtained in this wo rk fblack lines) w i th the 



Seifahrt et al 



(120071 ) (red lines) and by 



McElwain et al 



(120071 ) 



(blue lines). T he NIFS spectra were convolved by a Gaussian kernel to get similar spectral 



resolution as in 



Seifahrt et al. 



McElwain et al. 



(120071 ). namely 2500, 4000 and 4000 in JHK respectively. The 
(j2007l ) spectra are displayed at their full resolution of ~ 2000 in J and H. 
The main absorption/emission features are shown for each band. The spectra are normalized 
at 1.245, 1.68 and 2.29 /im respectively in JHK. 
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Fig. 6. — Comparison of GQ Lup b spectra with other known young and field brown dwarfs. 
The J-, H- and if-band spectra were respectively normalized with respect to their fiux at 
1.29, 1.68 and 2.29 ^m. The spectra were convolved with a Gaussian kernel to produce 
spectra at i? ^ 850. A constant was added to the different spectra for clarity. 
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Fig. 7. — Comparison of the GQ Lup b J-band spectra with spectra generated with the 
GAIA model at i? ^ 5000. Good fits found by a minimisation algorithm were obtained 
for Teff = 2300 - 2500 K and log g = 3.5 - 4.5 dex with a best fit for T^s = 2400 K and log g 
= 4.0 dex. The model spectra corresponding to the three gravity grid points at 2400 K are 
displayed as orange lines with the best fit value as the red line. The observed spectrum is 
depicted as the black lines. 
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Fig. 8. — Comparison of the GAIA synthetic spectra (red hnes) at Tcs = 2400 K and 
log g = 4.0 dex with GQ Lup b (black line) spectra in the H and K band at ~ 5000. 
The spectra were normalized at 1.68 and 2.29 /im. 
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Fig. 9. — Wide ba nd GAIA synthetic sp ectra compared t o pho tometr ic measurements 



taken in this work, in 



Marois et al 



fl2007f ). in 



McElwain et al 



(120071 ) and in 



Neuhauser et al. 



(120081 ). The model spectrum integrated flux in a given band is plot ted as a black dot. A s the 



measurements from this work and the CH4 and the V values from lMarois et al. 



mm were 



made by differential photometry wit 
variability of GQ Lup A reported by 



1 the host s t ar, th e error bars were set to include the 



Broeg et al. 



(120071 ). These error bars only reflect short 



term variations while long term fluctuations can also contribute to the discrepancy between 
predicted and measured values due to host star activity change. 
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Fig. 10. — GQ Lu p b position (diamond symbol) on the evolutionary tracks from 



Chabrier et al. 



(120001 ) that begins at 10^ years. The luminosity uncertainty is introduced 
by the lack of precision in the system distance measurement. The lower mass object of the 
eclipsing binary brown dwarf 2M0535-05 that has a dynamically determined mass of 36 ± 3 
Mjup is also plotted as a reference point. The evolutionary model underestimate the lumi- 
nosity of this object by a factor of ~ 1.5. Taking this systematic error into account leads to 



the blue error bars. 



